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Abstract. We propose the fundamental and two dimensional representa- 
tion of the Lorentz groups on a (3+l)-dimensional hypercubic lattice, from 
which representations of higher dimensions can be constructed. For the uni- 
tary representation of the discrete translation group we use the kernel of 
the Fourier transform. From the Dirac representation of the Lorentz group 
(including reflections) we derive in a natural way the wave equation on the 
lattice for spin 1/2 particles. Finally the induced representation of the dis- 
crete inhomogeneous Lorentz group is constructed by standard methods and 
its connection with the continuous case is discussed. 
Oh! 
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5h 1 1. Introduction 

a i 

The hypothesis of a discrete space and time has attracted the attention of physicists for 
different reasons [1] [3]: 

i) as a mathematical tool in order to get rid of the infinities in quantum field theories 
with the help of a cut-off in momentum space or, what is equivalent, a lattice for the 
space-time coordinates [7]. This model is similar to solid state physics where quantum 
fields are defined on grid points of a periodic crystal. 

ii) as a more realistic interpretation of lattice gauge theories, in which the space- 
time variables are constrained to discrete values due to some underlying structur resulting 
out of relations among fundamental processes, as Penrose, Finkelstein, Weizsaecker have 
proposed. [3] 

In this paper we address ourselves to an important problem of symmetries in lattice 
theories. In particular we study the consequences of restricting the continuous space- 
time variables to a discrete Minkowski space for the translations, rotations and Lorentz 
transformations. We discuss the possibility of maintaining the representation theory for 
these groups on the lattice. 
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Our paper is based on the standard theory of induced representations of the 
Poincare group restricted to discrete variables [13]. In momentum space there are three 
ways to construct induced representations: Mackey, Wigner and covariant states [10], 
based on the existence of a closed subgroup of a Lie group. In our case the closed 
subgroup is the cubic group with respect to the Lorentz transformations on the lattice. 
We have the advantage that the representations of the rotation group in two and three 
dimensions remain irreducible when restricted to the cubic group. Therefore all the 
arguments for the discrete case can be taken unchanged from the continuous case. 

In a preliminary version [15] of this work we have stressed the connection of Klein- 
Gordon, Dirac and Proca equation in discrete/continuous momentum and discrete space 
via Fourier transform. 

In this paper we emphasize the representation theory of the discrete translation, 
Lorentz and Poincare group in such a way that the wave equaton for spin 1/2 particles 
on the lattice emerges in a natural way from the Direc representation of the Lorentz 
group. 

In section 2 we describe an algorithm to construct all integral transformations of the 
complete Lorentz group based on the generators of some Coxeter group, and calculate the 
2-dimensional representations of this group that can be generalized to higher dimensional 
irreducible representations. 

In section 3 we review two unitary irreducible representations of the discrete trans- 
lation group and the cyclic group and use these representations as the kernel of two 
Fourier transform on the lattice that have become very helpful throughout the litera- 
ture. 

In section 4 we review the Dirac representations of the Lorentz group including 
space inversion and construct the Dirac wave equation in momentum space as the pro- 
jection operators that reduce the covariant states of the representation to the irreducible 
components. 

In section 5 we apply the Fourier transforms of section 3 to the Dirac equation in 
momentum space and obtain a difference equation for the Dirac and Klein-Gordon fields 
on the lattice. 

In section 6 we construct the induced representation of the Poincare group on 
the lattice using the Mackey- Wigner approach and discuss the irreducibility and orbit 
conditions of this representation. 

2. Fundamental and spin representation of the Lorentz group on the lattice 

An integral Lorentz transformation belongs to GL(4,Z) and leaves invariant the bilinear 
form 

r2 2 2 2 2 /i\ 

dj — ^0 1 *^2 * \ / 

According to Coxeter [1, page 47] all integral Lorentz transformations (including 
reflections) are obtained by combining the operations of permuting the spatial coordi- 
nates Xi,x 2 ,x s and changing the signs of any of the coordinates x ,Xi,x 2 ,x 3 together 
with the operation of adding the quantity x — X\ —x 2 — x 3 to each of the four coordinates 
of a point. 
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These operations can be described geometrically by the Weyl reflections on the 
planes perpendicular to the vectors 



oil = e x - e 2 , a 2 = e 2 - e 3 , a 3 = e 3 , a 4 = - (e + e 1 + e 2 + e 3 ) 



where {eo, &i, e 2 , e 3 ,} is an orthonormal basis. 
In matrix form these reflections are 
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These reflections generate a Coxeter group, the Dynkin diagram of which is the 
following 

4 4 

o o o o 

a 1 a 2 a 3 a 4 

Fig.l: Dynkin diagram for the Coxeter group which generates all integral Lorentz 
t r ansf ormat ions . 

Kac [2] has proved that Si, S 2 , S 3 , generate all the integral Lorentz tranforma- 
tions that keep invariant the upper half of the light cone. Note that Si, S 2 , S 3 generate 
the full cubic or octahedral group. 

The generator S\ can be used to factorize any integral Lorentz transformation that 
belongs to the complete (orthocronus) Lorentz group. Let 

( a e / 9 \ 
b 



c 
V d 



I 



be an integral matrix of determinant A = ±1, satisfying 

L*gL = g , g = diag (1, -1, -1, -1) 

and also a > 1. From 

a 2 _ b 2 _ c 2 _ d 2 = l 



(2) 



it follows that only one of b, c, d can be zero. Suppose a > 1. Then we apply Si, S 2 , S3 
to L from the left until b, c, d become non-positive integers. To the resulting matrix we 
apply S4. We get 

fa' e> f g'\ 
b' 



U 



C 

\d! 



J 
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with a' = 2a + b + c + d. Obviously 

(a + b + c + d) (a - b - c - d) = 1 - 26c - 2bd - 2cd < , 
therefore a + b + c + d<0 or 

2a + b + c + d = a' < a. 
By iteration of the same algorithm we get 

a>a' >a" > . . . > a (fe) > 1 . 

The last inequality is a consequence of the fact that L and S4 belong to the complete 
Lorentz group. Following this process we get an integral matrix with = 1 which is a 
combination of S 1: S 2: S 3 , giving all the 48 elements of the cubic group on the lattice. 

Therefore a general integral Lorentz transformation of the complete Lorentz group 
L can be decomposed as 



L = P?P°P"S A . . . S 4 P?P'P$ S A { S?S?S1 } (3) 

i2 ' j4 * 126 *l 126 )all permutations V 1 

where P\ = S1S2S3S2S1, P 2 = S2S3S2, P3 = S3 are matrices which change sign of 6, c, d 
and a, (3, 7, 8, e, (, 77, 9, 1 . . . = 0, 1 

In the continuous case the boosts characterize the quotient of the Lorentz group 
with respect to the subgroup of rotations. These continuous boosts take the vector 
(1, 0, 0, 0), stable under rotations, to any point of the unit hyperboloid in the forward light 
cone. In the discrete case the vector (1,0,0,0) is stable under the full cubic group acting 
on the space coordinates. The coset representatives of the integral Lorentz group with 
respect to the cubic subgroup are obtained by taking in equation (3) the quotient with 
respect to all integral elements of the cubic group denoted by {SiS^Snl 

l J all permutations. 

These coset representatives allow us to construct an algorithm to obtain any vector of 
the unit hyperboloid in the forward light cone from the vector (1, 0, 0, 0). Let (a, 6, c, d) 
be any vector with integral components satisfying a 2 — b 2 — c 2 — d 2 = 1, (a > 1). First 
we apply to this vector the parity operators Pi, P 2 or P3 defined in equation (3), in such 
a way that the components 6, c, d become non-positive numbers. To the resulting vector 
we apply the operator £4 obtaining a new vector (a', 6', c', d'), a' 2 — b' 2 — c' 2 — d' 2 = 1, 
with a' = 2a + b + c+d < a. Following the same procedure we obtain new integral vectors 
satisfying a > a' > a" > . . . > > 1. In the last step we obtain the vector (1, 0, 0, 0) 
as required. Taking the product of the generators used in this algorithm in inverse order 
we get the coset representative that takes the vector (1,0,0,0) to the vector (a,b,c,d), 
and we call this coset representative the integral boost. 

A particular case of integral Lorentz transformations are the integral Lorentz trans- 
formations without rotations. These can be obtained with the help of Cayley parameters 
[3] making n = p = q = and m, r, s, t integers. We can have the following cases: 

i) m 2 — r 2 — s 2 — t 2 = 1 



I 2m 2 — 1 2mr 2ms 2mt \ 

2mr 2r 2 + 1 2rs 2rt 

2ms 2rs 2s 2 + 1 2st 

\ 2mt 2rt 2st 2t 2 + 1 / 



(4) 
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ii) m 2 - r 2 - s 2 - t 2 = 2 
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(5) 



(6) 



(7) 



The solutions of the diophantine equations in the four cases are obtained by the 
application of all the Coxeter reflections as given in equation (3) to the vector (1, 0, 0, 0) 
in case i), to the vector (2,1,1,0) in case ii), to the vector (0,1,0,0) in case iii) and 
to the vector (0, 1, 1, 0) and (1, 1, 1, 1) in case iv). See Kac [2 page 70]. The four cases 
correspond to the integral Lorentz transformations given by Schild [1, page 42] restricted 
to the pure Lorentz transformations without rotations. 

Equations (4-7) can be considered also as particular cases of integral boost, that 
take the vector (1,0,0,0) to the vector defined by the first column of each of the four 
matrices. In the continuous case M0ller [4] has given a general boost that takes the 
vector (1, 0, 0, 0) to any vector on the unit hyperboloid. If we want to have a general 
matrix that takes the vector (1,0,0,0) to any integral vector (M, R, S,T) on the unit 
hyperboloid this is achieved by the square root of the first matrix of equation (4): 
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1+ 32 
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This matrix is the restriction of the general Lorentz transformation given by M0ller 
[4] to the values 

M = ^R = J^,S = ^,T 



P3 



m c 



m c 



m c 



m c 

satisfying M 2 - R 2 - S 2 - T 2 = 1. 

In position space the space-time coordinates of the lattice are integer numbers. 
They transform under integral Lorentz transformations into integral coordinates. The 
same is true for the increments Ax, 



'li- 



ft 



In momentum space the components of the four-momentum are not integer numbers 
but they can be constructed with the help of integral coordinates, namely, 



cAt Ax 



V((cAt) 2 -(Af) 2 ) 1/2 ' ((cAt) 2 -(Af) 2 ) 1/2 ; 



(9) 



If Ax^ transform under integral Lorentz transformations as a four- vector, the p^ will 
transform also as a four-vector because the denominator in each component is Lorentz 
invariant. 

In general there is no constraint between the values Ax^. But if we impose the 
condition (cAt) 2 — (Ax) 2 = 1 then the four-momentum can be considered an integral 
four vector multiplied by m c, m c(M, R, S, T), satisfying (M 2 - R 2 - S 2 - T 2 ) = 1 with 
M = cAt, R = Ax u S = Ax 2 , T = Ax 3 . 

Using the homomorphism between the groups SO (3,1) and SL(2,C) we obtain the 
representation of integral Lorentz transformations as 2-dimensional complex matrices. 
From the knowledge of the Cayley parameters [5] for a general element of the proper 
Lorentz group, we read off the matrix elements of the corresponding a G SL (2, C) 

1 fm + t + i(n — \) —p + r + i(q + s) \ . . 

^/~A~ \P + r + i(q — s) m — t + %{n + A) J ' 

with A = det a = m 2 — r 2 — s 2 — t 2 + m 2 + p 2 + q 2 — A 2 , 
m\ = nt + ps + qr. 

For instance we can calculate the 2-dimensional representation of the Coxeter re- 
flections Si, multiplied by the parity operator P (in order to get an element of the proper 
Lorentz group) identifying its matrix elements with the Lorentz matrix given in terms 
of Cayley parameters. Easy calculations give the unique solutions: 



a (PSi) 


1 


a (PS 2 ) 


1 


a (PS 3 ) 




a (PS A ) 





-1 -% 





% 1 
-1 -% 



l-i 
-l-i 2i 



(11) 
(12) 

(13) 

(14) 



From these matrices one can prove the expresions for the representation of the 
algebra for the Coxeter reflections, namely, 

a(SiSj) = a(PSi)a(PSj) , i^j , i,j = 1,2,3 

a (PSi) a (PSi) = -1 , i = 1,2,3 
a(PSi)a(PS 4 ) = a(SiS 4 ) , i = 1, 2, 3 (15) 
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The elements a (PSi) , i = 1,2,3 generate the so called octahedral binary or 
double group [14] excluding the parity. Together with a^PS^) they generate part of a 
binary version of the Coxeter group. 

The integral Lorentz transformations without rotations as given in eqs. (4) to (7) 
have two dimensional representations taking n = p = q = A = in eq. (10) and the 
choices 



i) m 2 — r 2 

ii) m 2 — r 2 

iii) m 2 — r 2 

iv) m 2 — r 2 



s 2 


-t 2 


= 1, 


s 2 


-t 2 


= 2, 


s 2 


-t 2 


= -1, 


s 2 


-t 2 


= -2. 



in eq. (4-7). In order to complete the picture we have to add the 2-dimensional repre- 
sentation of the matrix \f~L given in eq. (8) which turns out to be 

( fj\_ 1 ( M + l+T R-iS \ + 

a ^ =K = ^WTT){ R + ts m + i-t)=k- < 16 ) 

The 2x2 matrix representation of the discrete vector (M, R, S, T) , M 2 - R 2 - S 2 - 
T 2 = 1 can be obtained from the momentum in the rest system (1, 0, 0, 0) transformed 
by the matrix (16), namely, 

J 1 °)k + -( M + T R ~ iS ) (17) 

1 o i r ~( R+ts m-t ■ [n) 



3. Representations of the translation group on the lattice and Fourier trans- 
forms 

In the continuous case the unitary irreducible representations of the translation group 
are one dimensional and are given by the function 

D k (x) = exp {i2nxk) , x G R, (18) 

where A; is a continuous parameter that characterizes the representation. 

If we restrict the translations to a discrete set of points x = je, j G Z, the repre- 
sentation becomes 

D k (j) = exp (i2irkje) (19) 
where k is still continuous. This representation satisfies orthogonality relations 

f 1/2£ D k (j) D k (j') dk = -5 3f (20) 

and completeness relations 
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oo 



F(j) = [ 1/2£ exp (-2mkje) F (k) dk (23) 

J-l/2e 



E D k (j)D k 'U) = -S(k-k>). (21) 

j = -oo 

The parameter of the discrete translation group k is defined in the fundamental 

in 'Ye < k < Ye' 

From this representation we construct the Type I Fourier transform on the lattice 

oo 

F{k)= E exp (27rikje)F(j) (22) 

j = -oo 
»l/2e 
-l/2e 

If we generalize the discrete translation group to a (3+l)-dimensional cubic lattice 
T : {x II = j^e, fi = 0,l, 2 , 3, 3 ix e Z} 

with scalar product 

fc • x = fc*%e = (fc • j)e , (24) 
then the representations become 

D k (j li )=exp2m(k-j)e. (25) 

These representations satisfy orthogonal relations, completeness relations and 
Fourier transform analogous to the 1-dimensional case. 

We introduce the reciprocal group of discrete translations on the reciprocal lattice 
T^by 

The representations of the translation group in the (3 + l)-cubic lattice are invariant 
under the reciprocal group 

D k (j,) = D k + bR (j,) , (26) 

therefore, the parameters of the representations k are restricted to the fundamental 
domain 

The same property applies to the Fourier transform in (3+l)-dimensions 

F(k) = F (k + b R ) , b R e T R . (27) 
The irreducible representations of the reciprocal group can be written 

(b R ) = exp (2 7T if ■ b R ) , - G /2 < ^ <G /2 (28) 

satisfying 

f 72 (b R ) (b R ') d\ = e\n bR , (29) 



E (b R ) (b R ) = e*5 - O (30) 



Now we combine the irreducible representions of the discrete translations group in 
(3+l)-dimensions and the integral Lorentz transformations. 
Given a periodic function on the /c-space 

F (k) = F (k + b R ) , (31) 

it can be written in terms of discrete waves 

F(k) = ex P (2ttz (k ■ j) e) F . (32) 
U e r 

From the action of the group we have 

U A F(k) =P(A- 1 k) . 

Then 



U\F (k + b R )=F [A^k + A" 1 ^) 
But from the properties of the scalar product and using (32) we have 

exp (2ni (A-V • j) e) = exp 2m (b R ■ (A" 1 )^ e = exp 2ni (b R (A -1 )^/ 
for any integral Lorentz transformation. Finally, 

U A F (k + b R )=F (A^k) = U\F (k) . (33) 

Therefore, the periodicity of some function with respect to the reciprocal group is 
conserved under integral Lorentz transformations. 

Now we consider the one dimensional unitary irreducible representation of the cyclic 
group in one dimension of order A^: 

D m (j) = exp^mj , m,j = 0,l,---,N-l, (34) 

where j represents the space variable and m the label of the representation. 

This representation satisfies periodic boundary conditions with respect to the j 
variable 

D m (j + N) = D m (j), 
and also with respect to the m label 

D m+N (j) = D m (j) , 
and satisfies orthogonality relations 

1 £ D m (j) D m {f) = 5 mm , (35) 

JV m = 



and completeness relations 



-L £ D m (j) D m ' (j) = % (36) 
Iy£ j = o £ 



From this representation we construct a Type II Fourier transform on the lattice 

1 / 2-7T \ 

F (™) = ^ E ex P (^ mJ j F {j) (37) 

for any periodic function of discrete variable F(j + N) = F(j). This is called in the 
literature the Finite Fourier transform. 

The representation of the cyclic group in the one-dimensional lattice and the cor- 
responding Fourier transform can be generalized to the (3+l)-dimensional cubic lattice, 

2n 

D m (j») = expi—(m-j), (38) 

with (m • j) = m^j^ the bilinear form invariant under integral Lorentz transformations. 

Because of the periodic boundary conditions the label of the representation of the 
cyclic group on the lattice is constrained to the fundamental domain, namely, 

D m = D m (£„ + n,N) = D m , 

where < £ M < iV , ^eZ , G Z 
Symbolically 

D m (j) = D m (Z + f) = D m (Z), (39) 

where £ belongs to the fundamental domain and j c is any vector in the cubic lattice 
whose components are multiples of N, j c = {n^N}. The boundary conditions (39) are 
invariant under integral Lorentz transformations. From the definition of the group action 

U A f (x) = f (A" 1 *) , 

we have for a periodic function, periodic with respect to the space variable, 

(Uxf) (£ + f) = f (A^ + A- 1 /) = / (A-^) = (U A f) (0 , (40) 
where we have used the property 




due to the integral character of the Lorentz transformations. 

4. Dirac representation of the Lorentz group and covariant states 

Let L (a) be an element of the proper Lorentz group corresponding to the element a G 
SL (2, C) and I s the parity operator. One writes the components of four-momentum as 
a 2 x 2 matrix 

p = p fl o- f ,=p°a +p i o- l , (41) 
where <t = 1 and cr, are the Pauli matrices. 
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The transformations of p under parity and SL(2, C) are 

I s :p^p s = p°a - p l Oi = (det p) (py 1 , 

a : p — > apa + , (42) 
f -> (a + )"Va _1 . 

It follows 

/.L(a)/,- 1 = L ((Q+)" 1 ) . (43) 

Therefore the matrix (a + ) _1 gives another 2-dimensional representation of the 
Lorentz group non-equivalent to a G SL(2,C). In order to enlarge the proper Lorentz 
group with space reflection we take both representations a and (« + ) 1 . 

Let 7r = {I, I s } the space reflection group and a G SL(2,C), then the semidirect 
product 

SL(2,C) ® 7T 

with the multiplication law 



(a, 7r) (a', 7r') = irir') ifn — I, 

(a, it) (a', it') = (a(a' + ) if n — I s , 



(44) 
(45) 



form a group. 

This group has a 4-dimensional representation, particular elements of which are 



D(a,I) 



a 
(a 



+ i , D(eJ s 



a 
a 



(46) 



that satisfy 



D(e,I s )D(a,I)D(e,l; 1 )= D ((« + ) \ /) . (47) 

In this representation we could now construct D(Si) for the generators of the Cox- 
eter group. Using formulas (11)-(14) we get the 4-dimensional matrices: 



D(Si) 
D(S 4 ) 



a(PSi) N 

a(PSi) , 

[a+iPS^r 1 

a(PS 4 ) 



% = 1,2,3 



D(Si) D(Si) 



1 



i = 1,2,3,4 



DiSJ D(S 2 ) =±1, £>(S 2 ) L>(S 3 ) =±1, L>(S 3 ) L>(S 4 ) 



±1 
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With respect to this representation, a 4-component spinor ijj (p) in momentum space 
transforms as follows: 

U(a,I)^(p) = S(a,/)?(L- 1 (a)p), (48) 
U(e,I 8 )4(p) = D(e,I s )^(I s p). (49) 

Using a similarity transformation we obtain an equivalent representation 

D(a,ir) = MD{a,7r) M _1 



with 



In this representation 



M= 1 ( *o ^ 
y/2 \ -o"o o"o 



For this representation we can derive from eq. (46) the unitarity relation 

D+ (a, tt) D (e, I s ) D {a, n) = D (e, I,) . 
The new four-spinor 

transforms as 



U{a,I)i>(p) = D{a,I)il>(L- x {a)p), (52) 
U(e,I 8 )rl>(p) = D(e,I s )^(I s ,p), (53) 

and has an invariant scalar product due to the unitary relation given above. We call this 
representation the Dirac representation. 

The Dirac wave equation can be considered as a consequence of the relativistic 
invariance and irreducibility[10]: Under the restriction to SU(2), the first and second 
pair of components of the four-spinor transform according to spin 1/2. Irreducibility 
requires that one of these pairs should be eliminated. In the present discrete case we must 
replace the continuous group SU (2) by the binary octahedral group [14]. Fortunately the 
restriction of SU(2) to this discrete subgroup is irreducible [11]. This allows us to follow 
the steps of the continuous analysis [10]. In the rest system we want a projection operator 
that selects one irreducible representation of SU (2) out of the Dirac representation. This 
is achieved in the rest system by the projection operator: 
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In order to get the projection operator in an arbitrary system we apply D(k,I) 

given by eqs.(50) and (16). In the last equation if we identify M = R = S = 
-?L t=-&- , then, 

mo c rriQC ' ' 

Q^Q(k)= D- 1 (k, /) QD (k, I) = X - (I + W («)) , (55) 

where _ 

W (k) — - \ ^_ ) , ;_i , ; , ;_i (56) 

V ; 2\k + k-(k + k) \ -k + k-(k + k) 1 J K ' 



Using the identities 

(k + k) = —o^p^ > 



m c 

+ 1 7 

m c 



(57) 



we find 

W(k) = 7 %, (58) 

m c 

where 7 M are Dirac matrices with the realization 

7 ~ o -a ' 7 - 



7° = 7o , 7 J = ~ly 
Collecting these result we obtain the Dirac equation in momentum space 

Q{k)^{p) = 1 -{I + W{k))^{p) = ^{ P ) 

or 



(7%-m c/)V(p) = 0. (59) 

(An equivalent method can be used applying to the projection operator the Foldy- 
Wouthuysen transformation [9]) 

We apply the operator 7^^ + m c from the left to (59) and obtain the mass-shell 
condition 

(fp„ - m 2 c 2 ) 4> (p) = 0. (60) 

The Dirac equation is invariant under the group with elements (a, tt) defined before. 
In other words, if ip (p) is a solution of the Dirac equation, so is U (a, ir) ip (p). Put n = I. 
Then 
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U{a,I)^{p) = D{a,I)^{L- 1 p), 



D' 1 (a) Q (k) D (a) = Q (m) , 

Q H V (£~V) = V> (£~V) • 

Therefore 

Q( K )U(a,I)il)(p) = Q(K)D(a,I)*lj(L- 1 (a)p)=D(a,I)Q(Ka)ij(L- 1 (a)p) = 
= D(a,I)4>(L- 1 (a)p)=U(a,I)^(p) 

as required. For the space reflection 

Q (/,«) i) (i s p) = i) (i s p) , 

D- 1 (e,I s )Q(K)D(e,I s ) = Q(I sP ), 
U(e,I s )^(p) = D(e,I s )4>(I s p), 

we get 

Q( K )U(e,I s )if>(p) = Q( K )D(e,I s ),f>(I s p) = D(e,I a )Q(I a K)if>(I s p) = 
= D(e,I s )^(I s p) = U(e,I s )^(p) 

as required. Notice that, due to the relation between SU(2) and the binary octahedral 
group mentioned above, all properties of Dirac representation in continuous momentum 
space carry over to the discrete momentum space without modification. 

5. Dirac and Klein-Gordon equation on the lattice 

From the Dirac equation in momentum space (59) we can construct the wave equation 
in position space with the help of the Fourier transform we have introduced in section 3. 
We define the following difference operators 

A/(j)=/0' + l)-/0% A/(j) = i{/(j + l) + /(j)}, (61) 

v/0') = /0')-/0'-i), v/tf) = ^{/0') + /0'-i)}, (62) 

and the partial difference operators with respect to a function of several discrete variables 

A„/ (.7 M ) = / (j„ + <V) - / , 

Kf{j») = \{fU» + ^v) + fti»)}i (63) 
and similarly V„/ and VJ (j M ). 
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From these operators we construct 

k = \^ n ^ > = Jv, n v„ , (64) 

3 3 

^ + = n k > f = n^. (65) 

/x=0 /i=0 

From the Fourier transform we can derive the wave equation in lattice space. 
Type I Fourier transform: 

The kernel of the transform satisfies: 

1 2 

-A^ exp (27ri (k.j) e) — i- tan (irk^e) exp {2m{k.j)e) (66) 

We could apply the Fourier transform to the Dirac equation in momentum space 
(59) and would obtain the discrete wave equation. Instead we postulate a difference 
equation that in the limit goes to the continuous differential equation, namely, 

(^J-m c77 + )VO M ) = 0. (67) 

The kernel of the Fourier transform (25) or "plane wave" is a particular solution of 
(67) if it satisfies 

^7 M - tan rrk^e — m cj rj + exp 2m (k ■ j) e = 0. (68) 



Applying the operator tannk^e + m c from the left to the last equation we 



obtain 

4 

— (tan irk^e) (tan irk^e) — m\(? = 0, (69) 

which is the integrability condition for the solution of the wave equation. 

Now we multiply eq. (68) by some arbitrary (periodic) function ip (j M ) of discrete 
variable and sum for all j M . 

We get 

^7 M - tan nk^e — m cj ^ (k^) = 0, (70) 



where 

N-l 

4> (M = E i> U«) ^ ex P (^Trfc^e) (71) 

Jm=0 

is the Fourier transform of ip (j M ). 

If we compare eq. (70) with eq. (59) both are identical if we restrict the momentum 
p^ to the discrete values 

Pf, = -tanvr/c^e , k^ = -K- In - + (72) 
e 2me 1 - ^lep^ 

Instead of postulating Dirac equation on the lattice (67) we can deduce it in a 
natural way from the projection operator of the Dirac representation (59): first identify 
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the momentum p^ in this expression with the new variable as given in (72), then 
apply the Fourier transform of type I (23) and finally, using (66), obtain (67). 

Applying to the wave equation the operator i'"f fl 5~ +moci]~ from the left we obtain 
the discrete version of the Klein-Gordon equation in the lattice space 

(^--mgc 2 J7V)^0'A») = ° ( 73 ) 

a particular solution of which is again the "plane wave" (25) provided the integrability 

conditions (69) is satisfied 

Type II Fourier transform: 

The kernel of this Fourier transform satisfies: 

1 2ni .2 / 7r \ r 2ni , . . , . 

-\ exp — (j ■m)=i- tan I — m J A M exp — (j ■ m) (74) 

We postulate a difference equation on the lattice space the limit of which goes to 
the continuous wave equation, 

~ m o c V + ) i> tiv) = (75) 

The kernel of the Type II Fourier transforms (38) is a particular solution or "plane 
wave" of this equation (75) if it satisfies: 

/ 2 7r \ 2ni 

7 M - tan — m M - m c) r] + exp — (j ■ m) = (76) 
Applying the operator 7 M | tan j^m^ + m c from the left we get 

(tan (tan Jj m ») ~ m l° 2 = °> ( 77 ) 

which is the integrability condition for the "plane wave solution" of the wave equation 
on the lattice. 

Now we multiply eq. (71) by some arbitrary (periodic) function ip (j M ) of discrete 
variable and sum for all j M . We get 

(t M - tan - m c) $ {k„) = 0, (78) 

where 

N 2tt 

i> (K) = E i> tin) ^ +ex P i-j^^Jn ( 79 ) 

3n=0 

is the Fourier transform of ip (j M ). 

Both eq. (59) and eq. (78) are identical if we restrict the momentum p^ to the 
discrete values 

2 7r N 1 + -ispn 

Pn = - tan -m M , m M = — In \-. m M = 0, 1, • • • , N - 1 (80) 

e N 2ni 1 — pep^ 

From eq. (75) we could derive the Klein-Gordon equation on the lattice space as 
in eq. (73) with the integrability conditions (69) as before. 

As in the case of Fourier transform of type I we can deduce again the wave equation 
on the lattice in a natural manner from the projection operator (59): identify there the 
momentum p^ with the new variable given by (80), use the Fourier transform (37) of 
type II, together with (74) and obtain (75). 
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6. Induced representations of the discrete Poincare groups 

Let V\ = T^XgSO (3, 1) be the Poincare group restricted to the integral Lorentz trans- 
formations and discrete traslations on the lattice with the group composition 

(a, A) (a', A') = (a + Aa', AA') . (81) 

In order to construct irreducible representations we follow the standard method [13]: 

o 

(1) Choose an UIR, D k (a) , of the translation group T 4 

(2) Define the little group H e 5*0(3, 1) by the stability condition 

heH : D° k (h^aj = D° k (a) (82) 

(3) Choose an UIR D a of the little group H and construct for the group T 4 x s H 
the UIR 

D° k ' a (a, h) = D° k (a) <g> D a (h) (83) 

(4) Choose coset generators c of T 4 x s H in V\ constructed from the group action 

(a, A) c = c (a, h) (84) 

(5) Then the induced representations is: 

D% (a, A) = D k ' a (a, h) 5 ((c')" 1 (a, A) c, (a, h)) (85) 
and this is an UIR of V\. 

o 

Recall the construction of the massive UIR for the continuous case: If k is any 
vector inside the forward light cone, one can shift it on the orbit by a continuous Lorentz 

o 

transformation to the form k— m c (1,0,0,0), with the stability group H = 50(3). In 
the 2x2 matrix form this stability group becomes 50(2). 

o 

For the discrete Poincare group, we may choose k from the intersection of the 
Brillouin zone with the forward light cone. We can shift it on the orbit only by discrete 

o 

Lorentz transformations. Certainly we can choose k— m c (1, 0, 0, 0) within the Brillouin 
zone. Then the discrete stability group is the cubic group, and in the 2x2 matrix form 
it becomes the binary cubic group. In both cases the representations subduced from the 
continuous to the discrete little groups remain irreducible. 

For the coset representative c = (0, A) we can choose the integral Lorentz transfor- 

o 

mations A = L(k) that take k into an arbitrary integral vector of the unit hyperboloid. 
These transformations were defined in section 2 as integral boosts. The Dirac delta 
function in (85) is zero unless 

(a, h) = (0, L~ l (k')) (a, A) (0, L (k)) = (l- 1 (k') a, L~ l (k') AL (Jfe)) 

Substituting in eq. (85) with c, c' — > k, k' and using (82) we get 

Dl? k {a, A) = D° k (L- 1 (k 1 ) a) D a (l- 1 (k') AL (k)) (86) 

The spinor representation of the second factor is given with respect to the element 
L~ x (k') AL{k) that belongs to the little group, 577(2), or 50(3, R) respectively. These 
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representations of SU(2) corresponding to spin 1/2 or spin 1 state irreducible when 
restricted to the (binary) cubic group [11]. The first factor can be written: 

D° k (L- 1 (£/) a) = D k ' (a) (87) 



where k' = (L~ l (k')) T k are all the points that defined the UIR of the translation group 
and belong to the orbit on the dual translation group. This orbit is discrete in our 
analysis. 

We apply the analysis given in [13] for the semidirect product of the discrete trans- 
lation group and a point group on the lattice. The dual group of the translation group 
is given by all the points from the Brillouin zone. We wish to characterize the discrete 
orbit of the point group by function on the dual space. We formulate five conditions for 
these contraints: 

(1) they should vanish on the orbit points, (2) they should admit a periodic exten- 
sion on the /c-space, (3) the constraints must be Lorentz invariant, (4) the constraints 
should vanish only on the points of the orbit, (5) when the lattice spacing goes to zero, 
the difference equations in position space should go the continuous wave equation in 
Minkowski space. 

If we require only condition (l)-(4) the polynomials (k^k^ — k 2 ^ = vanish on 

and only on the discrete points of the orbit and the irreducible representation could be 
characterized by the functions (see Ref. [15], p. 192) 

f* (k) : (k»k, - k 2 ^ f* (k) = 0. (88) 

Nevertheless the difference equation (73) we have constructed in section 5, with the 
kernel of Fourier transform of Type I and II, does not lead to this constraint. We take 
a new approach. We choose the constraints we have derived in the Dirac representation 
in momentum space eq. (60) for the continuous case: 

(p% - m 2 c 2 ) i> (p) = (89) 

If we identify p M = | tan nk^e, with k in the Brillomin zone and e the lattice 
spacing and use the Fourier transform of Type I we obtain the difference equation (73) 

If we identify p M = - tan j^m^ in (89) and use the Fourier transform of Type II 
we obtain again the difference equation (73), the continuous limit of which leads to the 
continuous Klein-Gordon equation. 

The constraints (69) are periodic with respect to the /c M -space in the Type I and 
the constraints (77) are periodic with respect the m M -space in the Type II. 

Nevertheless when integral Lorentz transformations are applied to the components 
of the k^ or m M variables the new do not satisfy the constraint equations, and at the 
same time the constraints vanish at points not on the orbit. Therefore conditions (3) 
and (4) are violated, although they can be recovered in the asymptotic limit when e — > 

( - fc 2 U (M = Type I (90) 
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or in the limit e — > , Ne — > 2ir. 

(m»m^ - 1) tj) (m M ) = Type II. (91) 

The corresponding representation becomes irreducible and invariant. 

In both cases the wave equation in lattice space gives in the assymtotic limit the 
continuous Dirac equation. 
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